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We investigate, analytically and numerically, the fermion determinant of a new action on a 
(l+l)-dimensional Euclidean lattice. In this formulation the discrete chiral symmetry is preserved 
and the number of fermion components is a half of that of Kogut-Susskind. In particular, we show 
that our fermion determinant is real and positive for U(l) gauge group under specific conditions, 
which correspond to gauge conditions on the infinite lattice. It is also shown that the determinant 
is real and positive for SU(N) gauge group without any condition. 
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I. INTRODUCTION 



> 



As the Nielsen-Ninomiya theorem [Q states, we necessarily meet the difficulty of so-called fermion doubling problem 
when we formulate fermion fields on a lattice. In practical calculations, Wilson fermions have been widely used, 
where an additional term which vanishes in the naive continuum limit is introduced at the expense of the chiral 
\^ • symmetry. An alternative scheme was proposed by Kogut and Susskind In this scheme the chiral symmetry is 
maintained as discrete one and doubler fermions are regarded as fermions in other species. In (1 + D) dimensions, 
the Kogut-Susskind (KS) formalism describes a theory with 2^5~ degenerate quark flavors (2^+^ components). 
, Recently it has been shown that lattice fermionic actions satisfying the Ginsparg- Wilson relation [Q may provide a 
solution of the chirality problem In these attempts the modified chiral symmetry operator is used in stead of 75. 
The actions are local in the sense that the fermionic matrix are bounded by Ce^''''^'. However, it has been proved 
that actions with the Ginsparg- Wilson relation cannot be "ultralocal" From the practical point of view, the 
ultralocality (the couplings drop to zero beyond a finite number of lattice spacings) is also important. Thus ultralocal 
. fermionic actions with better features than, for example, the KS action is awaited, though not a final solution of the 
' chirality problem. 

, In the recent papers |@,||, we proposed a new type of fermionic action on a (1 -|- Z3)-dimensional lattice. The action 
is ultralocal and constructed so that fermion fields satisfy the bosonic type of dispersion relation. In this sense there 
are no extra poles in the propagator. We found that the minimal number of fermion components, dimensions of the 
spinor space, is 2^^^ in the Minkowski case and 2^ in the Euclidean case, which should be compared with 2^+-'^ of 
(jj the KS fermion. Furthermore our action has the discrete chiral symmetry as well. 

. It is much of interest to investigate the numerical feasibility of our new fermionic action. When dynamical fermions 
^ ' are included, the property of fermion determinants is crucial in numerical calculations. For example, some methods 
• i-H : proposed to treat fermionic freedoms rely on reality or positivity of the fermion determinants |^ . 

In this paper we report the analytical and numerical results on the fermion determinants of our new action in (1 + 1) 
H [ dimensions. Our main concern is on U(l) gauge group, but some results on SU(N) gauge group are also presented. 
- - '■ In the case of U(l) gauge group, we will see that calculations with specific conditions for temporal link variables are 
stable and satisfactory though results without the conditions are unstable. The reason why we need those conditions 
will be discussed in detail. 

In Sec. 2, we recapitulate our formalism for later convenience. The analytical and numerical results on fermion 
determinants will be presented in Sec. 3, which is followed by the summary. 
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II. NEW FERMIONIC ACTION 



In the previous paper we proposed a new fermionic action on the Euclidean lattice. Though the action respects 
the discrete chiral symmetry like one in the KS action, fermion fields in this action have 2^ components in (1 + D) 
dimensions, which should be compared with 2^+^ in the case of KS fermions. In this section we briefiy sketch our 
formalism for later convenience. 

The action can be written with the fermion matrix A as 

= ^ ■0m^m,nV'n, (2.1) 
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where the summation is over lattice points and spinor indices, and our fermion matrix is defined by 

A = 1 - sIUe. (2.2) 
Here Ue is the Euclidean time evolution operator and is the unit shift operator defined as 

Si,^P{x°,x\...,x>',...,x'') =ij{x°,x\...,x^ + l,...,x'') (m = 0,1,...,£)). (2.3) 
We required that the propagator has no extra poles and found that U e has the form 

D 



E = l-Y,'f {iXi[Si- Si) +{!-¥,) [s,-2 + s})}, (2.4) 

i=l 

where is the ratio of the temporal lattice constant to the spatial one and X's and Y's, which are matrices with 
respect to spinor indices, should satisfy the following algebra: 

{Xi , Xj } = — Sij , 

{Xi,Yj} = 0^ (2.5) 
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{Yi,Yj} = 2(^— % + l 

where i and j run from 1 to D. The matrix 2{5ij/rE + 1) is positive definite for any positive r^, therefore X's and 
F's can be assumed hermitian, 

xl=Xi, y}=Y,. (2.6) 

The matrices X's and F's are written in terms of the Clifford algebra Fi, . . . , V^u as 

p. ^ 

Xi = Y, = y^aijFD+j, (2.7) 

V' E 



where 

rj = r„ {r„r,} = 2<5,, (»,j = i,...,2D), (2.8) 

and aij are certain real constants. The dimension of the irreducible representation for F's is 2^ and accordingly tjj 
has 2^ components. 

Now we give some useful properties of the time evolution operator Ue- We can see immediately that 

4 = Ue, (2.9) 
as X's and F's are hermitian. Since the inverse of Ue is given by 

= 1 - E T {-'^' (^^ ^0 + + {Si-2 + Sj) } , (2.10) 

we find that Ue is related to its inverse as: 

T2D+iUeT2D+1 = Ue\ (2.11) 

where 

r2D+i = iri,r2 . . .r2D, (2.12) 

r2D+i=r2D+i, r2^+i = i, {ri,r2D+i} = o ii = i,---,2D). (2.13) 

The interaction of the fermion with gauge fields is introduced by replacing the unit shift operators by covariant 
ones: 



where /t is the unit vector alon g the /i'th direction, and Ux.y is a hn k va riable connecting sites x and y. 
The fermion matrix Eq.(2.2) and the time evolution operator Eq.(2.4) become 



A{x) = 1-SI{x)Ue{x), 



(2.14) 



(2.15) 



and 



Ue{x) = 1 - E Y {'^' (^'^"^^ ~ ^'^^"^0 + " (^^("^^ " 2 + ^' ("^0 } 



(2.16) 



With gauge fields coupled to the fermions, however, Eqs.(2.10) and (2.11) do not hold any more in arbitrary 
dimensions except for (1 + 1) dimensions. In (1 + 1) dimensions these relations are 



(2.17) 



and 



T:iUE{x)T^ = U^\x), 



(2.18) 



which are used in the next section. 



III. ANALYTICAL AND NUMERICAL RESULTS OF OUR FERMION DETERMINANT IN (1 + 1) 

DIMENSIONS 



A. Reality of the determinant 

In this section we show the analytical and numerical results of our fermion determinant in the (1 + l)-dimensional 
case. 

First, we investigate the determinant of the unit shift operator 6*0(2;) in (1 -\- D) dimensions. Let us denote the link 
variables by Ux.x+p, which live on the links connecting two neighboring lattice sites x and x + fi. Now we consider the 
case of the Abelian gauge group U(l), and they can be written in the form 



U. 



X.X-\-}.L 



(3.1) 



where Of^{x) is restricted to the compact domain [0, 2tt). 

Then, 150(2;) can be constructed as a NtNxNs x NtNxNs matrix where Nt and Nx is the size of the lattice in the 
Euclidean time and the spatial directions, respectively, and Ns is the number of spinor components. >S'o(2;) is diagonal 
in the spatial and spinor space and consists of block Nt x Nt matrices Q{x) belonging to the Euclidean time space. 
Thus we can write the determinant of £'0(2;) as follows: 



det Sq{x) 



Q(i) 



Q(2) 



(3.2) 



where Q{x) for each x appear Ns times as 5*0(2;) is unity in the spinor space. 
The determinant of Q{x) is easily calculated by noticing its matrix form: 
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det Q{x) 



U(^2,X),{3,X) 



'U(Nt,X),(1,X) 



U(^Nt-i,X),(Nt,X) 





(3.3) 



^U(^l^X),(2,X)U(^2,X),{3,X) ■ ■ ■ {~U(^Nt.X),{1,X)] 



(-)^-exp<^z^0o(t,a;) 



(3.4) 



We used the anti-periodic boundary condition for the time direction, which is represented by a negative sign in the 
lower-left corner of the matrix Q. 

As a result, in {\ + D) dimensions we find 



det5o(a:) =n 



X 



Nt 



(-)^-cxp *^0o(i,a^) 



i = l 



exp|iiV,^^^o(t,a;)| =e^", 



where the angular variable a is defined by 



(3.5) 



(3.6) 



t X 



From now we as sume D = 1 and show that the determinant of the Euclide an t ime evolution operator Ue{x) is 
unity. By Eq.( 2.1S| ) we easily find detJ7_B(a;) = ±1. Since X's and F's obey Eq.(2/7) and in the limit ^ 

teX, 0, teY, 0, (3.7) 

det Ue{x) is equal to 1 at r^; = 0. From the continuity of det Ue{x) with respect to te > 0, we conclude 

dctC/B(x) = 1, (3.8) 

in the (1 -f l)-dimensional case. 



Collecting the above results, the relations (2.18) and (3.8), and the hermiticity of Ue{x), we can obtain a relation 
on the phase of our fermion determinant in (1 + 1) dimensions as follows: 

(detA(x))* = det (l - Sl^xjUEix)^ = det (1 - Ue{x)Sq{x)) 



det 



Ue{x) det So{x) det {sl{x)UE{xy^ 



(_)dim. of Adet5o(x)detA(a;). 



(3.9) 



Since the number of components of our fermion is Ns — 2 and consequently the fermion matrix A(a;) always has an 
even number of dimensions, we get 



(detA(a;))* = e'"dctA(a;). 



(3.10) 



Accordingly in the case of (1 + 1) dimensions with U(l) gauge fields, our fermion determinant is real under the 
condition 



a = 27rn (n: integer), 



(3.11) 



where a is defined in Eq.( |3.6D . The condition specified by Eq.( |3.1l| ) will be called " GT-condition" in this paper, 
where the temporal link variables are globally constrained as X^x^ol^) = T^'n and on the infinite lattice the condition 
is achieved by a gauge transformation. We also define " T-condition" as 9o{x) = 0, which corresponds to the temporal 
gauge condition on the infinite lattice. 
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At a glance, it seems strange that the determinant has a-depcndence, because a can be gauged away. However, 
this is not true for a finite lattice, since all gauge transformations on the infinite lattice are not allowed on a finite 
lattice with periodic boundary conditions. Thus, the dependence comes from finiteness of the lattice. For example, 
J2t^o{t,x) is invariant under any gauge transformation on a finite lattice with periodic boundary conditions. The 
temporal gauge condition is certainly a consistent gauge in an infinite system, but not in a finite system. However, 
this restriction is just an artifact arising when we try to approximate an infinite system by a finite one. We could think 
in the following way. Given an infinite system, we impose the temporal gauge condition, which is legitimate. Then 
we take an finite sub-system to approximate the original total system. The difference is due to size and boundary 
effects and would disappear in the infinite volume limit. The same is also true for the GT-condition. It should be 
noted that the fermion determinant is invariant under gauge transformations on a finite lattice, but the determinant 
does depend on our "gauge conditions": the T- or the GT-condition. 
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FIG. 1. The distributions in the complex plane of our fermion determinants for each configuration of 2000 Monte Carlo 
iterations after getting good enough equilibrium, i.e. after 2000 iterations, on a 8 x 8 lattice at /? = 1.5. Plots (a),(b) and (c) 
correspond to under no conditions, the T- and the GT-condition, respectively. 
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Now, we will give the numerical results in the (1 + l)-dimensional U(l) gauge theory with our fermion action. In the 
following numerical simulations link variables are updated by the Metropolis method and determinants are calculated 
by the LU decomposition. So there are no systematic errors in the determinants. The way to generate a sequence of 
configurations in the Monte Carlo computation is as follows. In the T-condition, we fix Uo{x) — 1 everywhere. After 
this, all temporal link variables are not changed and we only update the link variables of spatial direction using the 
Metropolis method. In the GT-condition, we prepare a configuration satisfying '^^Oa{x) — 7m. Then a temporal 
link variable e*^°^^^ is replaced by e^^°^^^^^^, where % is a random number between — tt and tt. At the same time 
another temporal link variable, which is chosen randomly, is multiplied by e~*-^. And we use the above configuration 
as the trial one in a Metropolis acceptance test. It can be shown that this procedure satisfies the micro-reversibility 
requirement, and therefore also the detailed balance condition. 

In Fig|^(a) - |l|(c) we show the distribution of the fermion determinants in the complex plane for some (3. The 
statistics is 2000 thermalizations and 2000 measurements. In Fig.|l|(a), no conditions are imposed on temporal link 
variables. Fig|^(b) is t he re sult in the T-condition. In Fig.|l|(c), the GT-condition with n = is imposed, namely the 
condition a = in Eq.(3J) is always kept in Monte Carlo updates. The distribution in Fig|^(a) has a doughnut-like 
structure with a center around the origin, so that we can expect that the convergence for any observation is very 
poor. On the other hand, Fig||(b) and]l|(c) show that the determinants in the T- and the GT-condition are real as 
expected and also positive. The distribution in the imaginary direction is due to numerical errors (note the difference 
in the scales of real and imaginary parts). The positivity of determinants is important from the numerical point of 
view and will be discussed later. 

To test the convergence in the above three types of conditions, we measured the expectation value of a plaquette 
value as a function of /3 using the following formula: 



dUdijdij^'e- 



dUWp dctA{x)e~''^ 
dU det A(a;)e"'^'' 



{Wp detAjx))^ 
(detA(a:))o 



(3.12) 



Here Sg is the usual action for link fields, and dip, dip'^ and dU stand for Ylx a dipaix), Yix a ^i^ai^) ^^'^ Tlx dUx,x+fi, 
respectively. 





O O No conditions 
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FIG. 2. The real part of the averaged plaquette vs /3 on a 8 x 8 lattice. 
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In Fig.g we can see that the expectation values in the T- and the GT-condition display gentle curves, while one 
without any conditions is spiky especially for small (3. Here, the plotted points are the average over 10 results, each of 
which is obtained by 2000 Monte Carlo iterations at each /3. And the error bars are evaluated by using the standard 
deviation of the 10 results. If the error bars are not displayed, they are not visible at this scale. 

T he po orness of the convergence without any conditions can be traced back to the behavior of the denominator in 
Eq.( |3.12 ). The expectation value (det A(x))o not only takes complex values but also suddenly increases or decreases 
during sampling. Fig.H shows the absolute value of the averaged fermion determinant as a function of the update 
iterations for each condition. We find that the convergence in the case of no conditions is ill while in the cases of 
other two conditions very fine. Generally the more iterations are expected to improve the convergence. Without any 
conditions, however, this is not the case since the expectation value of the determinant must be in the empty hole of 
a doughnut-like structure in Fig|^(a) . Thus it is difficult to improve the convergence for any expectation value within 
a reasonable number of iterations. 

Now what is the difference between the T- and the GT-condition? We expect that the GT-condition is superior 
to the T-condition since the former condition is much weaker condition than the latter. To see this, we show in 
Figs.^(a) and ^(b) the averaged plaquette value in the pure U(l) gauge theory as a function of /3 for three types of 
conditions. We find the line without any conditions and one with the GT-condition are very close to each other, but 
the line calculated in the T-condition is slightly upper than those. As expected, this difference comes from the size 
effects in each condition and obviously tends to decrease as the lattice size becomes larger. Thus we conclude that 
the GT-condition has better feature on a finite lattice: the fermion determinant is positive and the finite-size effects 
are smaller. 
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FIG. 3. The absolute value of the averaged fermion determinant vs Monte Carlo iterations. The data were computed using 
the same configurations as in Fig. 1(a) - 1(c) 
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FIG. 4. The averaged plaquette vs /? on a (a) 4 x 4 and (b) 8x8 lattice. 



B. Spectrum of the fermion matrix 



Before discussing the positivity of the fermion determinant, we study the spectrum of A{x). First, we introduce a 
discrete rotational symmetry in the complex plane for eigenvalues of S'Q(a::)[/B(a;). Defining {T)t^t' — <2xp{i27Tt/NT)St.t' , 
we find 

TSI{x)Ue{x)T-^ = exp (^i^^ sI{x)Ue{x). (3.13) 

This relation implies that if A is some eigenvalue of SI{x)Ue{x), then exp (i27r/A^T) A is also its eigenvalue. Next 
using the following relation 



det {Ue{x)Sq{x) - A*l) = det {Ue{x)So{x)) det - X*Sl{x)UE{xy 

= det {Ue{x)So{x)) det (l - X* sI{x)Ue{x)^ , (3.14) 



we rewrite the eigenvalue equation 



into the form: 



det 



(^sI{x)Ue{x)~ \l) ^0 (3.15) 



det 



(^Sl{x)UEix)-\*-^l) =0. (3.16) 



Accordingly, from Eqs.( 3.1^ ) and ( |3.16| ) if A is some eigenvalue of sI{x)Ue{x), then A* ^ is also its eigenvalue. 

Let us look at the numerical results of the spectrum of our fermion matrix A(a;) = 1 — sI{x)Ue{x) in the (1 + 1)- 
dimensional U(l) theory. Figs.^(a) - ^(c) display the spectrum in the case of no conditions, the T- and GT-condition, 
respectively. From the figures we confirm the two properties of the spectrum discussed above for each condition. The 
spectrum in the GT-condition is so similar to the one with no conditions that one couldn't distinguish them at a 
glance. It is interesting that the determinant is real only in the former case. 

Now, using above two properties for Sq{x)Ue{x), we discuss the positivity of our fermion determinant. We will 
give a plausible reason for the positivity not a complete proof. 
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First we write the determinant as follows: 



det 



{l-Slix)UEix)) 



1 



(1-Ai)(l-Ai*)(l-A2)(l-A2*) 



(-) 



NtNx 



Ai' 



A. 



(3.17) 



In Eq.(3.17), the denominator must be real, since the numerator is positive and detA(a;) is real under the T- or 
the GT-condition as shown before. The denominator is a continuous function of the background configuration and 
furthermore it is never vanished, because sI{x)Ue{x) has always an inverse. Once the denominator is positive for 
some background configuration, it keeps on having a positive value for any configuration, if one configuration can be 
transformed continuously to the other within the T- or the GT-condition. 



(a) No conditions 
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(b) T-condition 
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(c) GT-condition 
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FIG. 5. The spectrum of A(a;) in the complex plane for a 8 x 8 lattice in the case of (a) no conditions, (b) the T- and (c) 
the GT-condition, respectively. 
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Next, we show that our fermion determinant is positive when all link variables are set unity. In this case sIUe can 



be easily diagonalized in the momentum space. The eigenvalues are expressed by the eigenvalues Xum Ue as 



2n + 1 

At/mexp(i— — tt) (n 0, . . . ,iVT - 1). (3.18) 



detA= [] (l-At/„,exp(z^^^7r) ) (3.19) 



Therefore 

/ / 2 

det A = II I 1 — Arr... exn I i- 

Nt 

is clearly positive, as Xijm can be shown positive. 

The configuration with all link variables unity clearly satisfies the T- and the GT-condition. This seems to complete 
our proof for the positivity of the determinant. It should be noticed that we might be faced with configurations where 



some A with |A| = 1 is not degenerated in spite of the symmetry A and A* . In this case the relation Eq.(3.17) does 
not hold and our proof fails. However, such configurations are very hard to happen. 



C. SU(N) gauge group 

When the gauge group is SU(N), we can also make simila r di scussion to the U(l) gauge group case and find the 
positivity of our fermion determinant. In this case, from Eq.(|3.3[) we immediately see 



(3.20) 



detQ(a;) = ± \U{i,x),{2,x) \ \U(2^x),{'i,x) \ ■ ■ ■ |f^(JVr-i,a;),(o,a;) | 
= ±1, 



so we find det 5*0(2;) = 1 without any conditions. Consequently, the determinant of our fermion matrix is always real 
in the case of SU(N) gauge fields. And we find that our fermion determinant is positive in the (1 + l)-dimensional 
SU(N), since the two symmetries of the spectrum of h.{x) discussed above are also satisfied. In fact we can confirm 
the symmetries from the spectra shown in Figs.||(a) and ^(b) for SU(2) and SU(3), respectively. And, in Figsj^(a) 
and|^(b), we display the numerical results for the distribution of our fermion determinant in the (1 + l)-dimensional 
SU(2) and SU(3) gauge theories. The results show that fermion determinants are real and positive in both cases. 
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FIG. 6. The spectrum of A(a;) in the complex plane for a 
gauge groups, respectively. 
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lattice. Plots (a) and (b) correspond to SU(2) and SU(3) 
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FIG. 7. The distributions of our fermion determinant in the complex plane for each configuration of 2000 Monte Carlo 
iterations on a 8 x 8 lattice ai j3 — 2.0. Plots (a) and (b) correspond to SU(2) and SU(3) gauge groups, respectively. 

We conclude that our fermion determinant det A(a;) in (l+l)-dimensions is real and positive for U(l) gauge group 
under the T- or the GT-condition and for SU(N) gauge group without any conditions, though we have only a plausible 
reason for the positivity. 

IV. DISCUSSION AND SUMMARY 

When we applied our new action without any conditions to U(l) gauge theory on a lattice, we were faced with 
the problem of convergence in Monte Carlo simulation. In this note, we showed that we could avoid this problem by 
imposing the T- or the GT-condition. In order to make this situation clear, as an example, we study the propagator 
of the fermi field. 



di}di}UUil}{y)i}\y')e-^o-^f 



(4.1) 



Integrating in Eq.(4.1) with respect to ^-nd ?/;^(a^) we have 



mv)i^\y')) 



j dC/det (l 


- sI{x)Ue{x)) 


(l - Slix)UEix)) 


-1 
v,v' 




f dU det (l - Sl{x)UEix) 





(4.2) 



We introduce new variables 8 and 9o{x) instead of link variables = e'^^o(x)^ 



^0(2:) = 6 + ^0(2:), 



(4.3) 



and using the following equation 



lld9o{x)d&sij2^"(''n 

X \ X / 

Y\deo{x)dQ5 {"^eoix) - NtNxO 

X \ X ) 
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we get 

J dUdQS i 9oix) J det (l - e'®5o(.T) ^^(a;)) (l - e*®S'o(a;) ^^(a;)) \ e"^" 



J dU deS (Y^^oixU dot (l - e'®S'o(x)V£;(x)) e~ 



(4.4) 



where the tilde represents replacing ^0(2;) by ^0(2;). 

The phase of determinant det(l — e^^ So{x)'' Ue{x)) takes any value between and 2tt as is shown in Fig.|l|(a). The 
summation of det(l — e^^" So{x)^ Ue{x)) o ver a sequence 0n(n =1,2,...) which is chosen at random are canceled out 



accidentally, then the denominator of Eq.(4.4) becomes very small. This is the origin of unstable behavior in Monte 
Carlo simulation. Indeed in Sec. 3, under the T- or the GT-condition which ensures that the variable Q is fixed zero, 
we proved the determinant det(l — Sq{x)Ue{x)) in the (1 + l)-dimensional lattice to be real for all configurations of 
link variables and to be positive for most ones. 

We have another reason for fixing O, namely, imposing some condition like the T- or the GT-condition. Since the 



integrand of the numerator in Eq.(4.4) is equal to the cofactor of matrix (1 — e So{x) Ue{x)) whose elements are 



linear functions of e , it is a polynomial of e . After we integrate it with respect to 0, all terms of this polynomial 



vanish besides constant terms. For the denominator of Eq.(4.4) we can say the same thing as the above, thus we have 

{i>{y)i^\y'))^ly,y'. (4.5) 

which is an undesired result. Contrarily, if we choose some condition, like the T- or the GT-condition, we may avoid 
this trouble. 

When we try applying this fermion to the SU(N) lattice gauge theory, we can expect that there is no necessity 
for imposing some condition because the element like e'® does not belong to the group SU(N). In fact, it is shown 
that the fermion determinant is real for all configurations and that it is positive for most configurations. And the 
numerical simulation for SU(2) and SU(3) gauge theory shows that the fermion determinants are real and positive. 
We will discuss the apphcation of our fermions to SU(N) lattice gauge theory in higher dimensions elsewhere. 
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